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We propose a new method to describe three-body breakups of nuclei, in which the 
Lippmann-Schwinger equation is solved combining with the complex scaling method. The 
complex-scaled solutions of the Lippmann-Schwinger equation (CSLS) enables us to treat 
boundary conditions of many-body open channels correctly and to describe a many-body 
breakup amplitude from the ground state. The Coulomb breakup cross section from the 
6 He ground state into 4 He+n+n three-body decaying states as a function of the total exci- 
tation energy is calculated by using CSLS, and the result well reproduces the experimental 
data. Furthermore, the two-dimensional energy distribution of the El transition strength 
. is obtained and an importance of the 5 He(3/2~) resonance is confirmed. It is shown that 

^ ■ CSLS is a promising method to investigate correlations of subsystems in three-body breakup 

reactions of the weakly-bound nuclei. 

in 

§1. Introduction 

in 

<^ | A neutron halo structure is one of the most interesting topics in physics of 

neutron-rich nuclei. In particular, the two-neutron halo structure observed in the 
Borromean systems such as 6 He and 11 Li, where any binary subsystem does not 
have a bound state, has attracted much attention and has been studied by many au- 
^ 1 thors. 1 ) -3 ) Theoretically, many works have been performed to understand a binding 

mechanism of these nuclei based on core+n+n three-body models, and an impor- 
tance of two-neutron correlations has been pointed out. 2 ) -7 ) Recently, it has been 
discussed how to clarify the internal correlations of core-ra and n-n subsystems in 
the two-neutron halo nuclei from observables. 8 )' 9 ) 

Coulomb breakup reactions using a high-Z target such as Pb have been con- 
sidered as useful tools to investigate the weakly-bound halo nuclei. For 6 He, the 
Coulomb breakup cross sections were measured by GSI 11 ) and MSU 12 ) groups. For 
Li, there were three sets of data measured at MSU, 13 ) RIKEN 14 ) and GSI, 15 ) and 
recently, a new measurement at RIKEN 16 ) was reported by Nakamura et al. Through 
those data, we can obtain the understanding not only of ground state properties, but 
also of breakup reaction mechanisms of halo nuclei. Especially, for two-neutron halo 
nuclei, the observed cross section is expected to give some important information on 
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internal correlations of core-n and n-n subsystems. To understand the correlations 
of subsystems, it is necessary to investigate the Coulomb breakup reaction based on 
a reliable theoretical approach. 

In our previous studies, we have successfully described the Coulomb breakup 
reactions of 6 He and 11 Li by using the extended core+n+n three-body model and 
the complex scaling method (CSM). 17 ^ 19 ) In these analyses, the cross sections were 
calculated using the response function method (RFM) combined with CSM, which 
is based on the linear response dominated by the El strength. For 6 He, the strength 
distribution is found to have a peak at around 1 MeV, which is dominated by the 
transition into the 5 He(3/2~)+n two-body continuum states. 17 ) This result indicates 
that the sequential breakup process via the 5 He(3/2~)+n components is important, 
and it is shown in Ref. 17) that the low energy peak originates from the threshold ef- 
fect reflecting the halo structure of the ground state. In the n Li breakup case, it was 
shown that the non-resonant three-body continuum states of 9 Li+n+n give a com- 
parable contribution with the sequential process via 10 Li+n in the cross section. 19 ) 
For both cases, our results well reproduced the observed breakup cross sections with 
respect to total excitation energies. 

In the previous analysis, correlations of subsystems were investigated by sep- 
arating the transition strength into the components of resonant and non-resonant 
continuum states. This separation of the strength is useful when we discuss the 
effects of resonant and non-resonant continuum components on the structure of the 
strength, and further clarify how much the total strength is exhausted by the indi- 
vidual strength. While the total strength obtained using RMF reproduces the ex- 
perimental observable, however, the separated strength does not correspond to the 
observable directly since the experimental data always contains both contributions 
of resonant and non-resonant continuum states. Then, in this study, we consider 
another approach, which describes the observables exhibiting the information on in- 
ternal correlations in a three-body breakup. To extract internal correlations in a 
three-body system from the observable, it is essential to describe the physical quan- 
tities as function of relative energies and momenta in binary subsystems. In fact, 
experimentally, the breakup cross section was reanalyzed as a function of subsystem 
energies to understand correlations in two-neutron halo nuclei. 8 ) 

Theoretically, it is a difficult problem to describe physical quantities of three- 
body breakups with composite particles having internal structures. The standard 
methods such as the Faddeev, of course, work well when we handle a simple three- 
body scattering with point particles. However, it is difficult to apply them to the 
composite particle case. Therefore, an alternative method is needed. 

For the two-body case, in Ref. 20), it is shown that we can calculate the exact 
scattering amplitude by using the formal solution of the Lippmann-Schwinger equa- 
tion (LS Eq.) with the complex-scaled Green's function even if we handle a scattering 
problem with composite particles. It is noticed that the Green's function in Ref. 20) 
is constructed by discretized eigenstates of the complex-scaled Hamiltonian, which 
are solved in the same manner as bound state cases, and satisfies correct boundary 
conditions without any explicit enforcement of boundary conditions. Furthermore, it 
is shown that this complex-scaled Green's function also works to describe the three- 
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body breakups. 17 ^ It indicates that we can easily apply the procedure in Ref. 20) to 
three-body cases. 

Additionally, the formal solution of the LS Eq. is useful to describe physical 
quantities of three-body breakups as functions of relative energies in subsystems, 
because it is represented by a solution of an asymptotic Hamiltonian, namely, a 
plane wave of a three-body system. 

The purpose of this work is to extend the theoretical approach in Ref. 20) to 
three-body systems and propose a new method which can evaluate scattering ob- 
servables as functions of energies in subsystems in three-body breakup reactions. 
In this paper, we apply this method to the Coulomb breakup reaction of 6 He and 
show that this method is capable of investigating internal correlations of subsys- 
tems in three-body decaying systems. The reliability of this method is shown by 
calculating the Coulomb breakup cross section of He. Furthermore, we evaluate the 
two-dimensional energy distributions of the El transition strength associated with 
the subsystems in 6 He, which is useful to investigate the internal correlations. In 
particular, the importance of the 5 He(3/2~) resonance in the final states is confirmed. 

This paper is organized as follows. In § 2, we give an explanation of our method 
to describe three-body scattering states as functions of subsystem energies. In § 3, 
we show the obtained results of the Coulomb breakup reaction of 6 He, and discuss 
the reliability of our method and the correlations of subsystems seen in this reaction. 
The last section, § 4, contains a brief summary. 

§2. Complex-scaled solution of Lippmann-Schwinger equation 
for three-body breakup 

In this section, we explain our new method to describe the three-body Coulomb 
breakup reaction in an energy representation of subsystems. Before describing our 
method, we give brief explanations of the He+n+n three-body model of 6 He and 
CSM in 2.1 and 2.2, respectively. In 2.3, we describe the formalism of our method 
named as the complex-scaled solutions of the Lippmann-Schwinger equation (CSLS). 

2.1. 4 He+n+n model of ' 6 He 

We first explain the 4 He+n+n three-body model of 6 He briefly. More detailed 
explanation is given in Ref. 17). In this model, we describe the 4 He core as the 
(0s) configuration, whose oscillator length b c is taken as 1.4 fm to reproduce the 
charge radius of 4 He. In order to analyze the breakup reactions, it is important to 
reproduce a threshold energy for each open channel and scattering properties of every 
subsystem correctly. Hence, we employ the orthogonality condition model (OCM), 21 ) 
in which we can use the reliable Hamiltonian whose inter-cluster potentials satisfy 
the conditions mentioned above. 

We solve the following OCM equation for the relative wave function x J7r of the 
4 He+n+n system; 



H X J \nn) = E X J \nn), 



(2-1) 
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where the Hamiltonian for the relative motion is expressed as 

3 2 

H = J2k-T G + J2 Van(n) + Vnn + t„ + ^PF- (2-2) 

i=l i=l 

The operators ti and Tq describe a kinetic energy of each cluster and a center-of- 
mass motion of a three-body system, respectively, and rj (i = 1 or 2) represents 
a relative coordinate between 4 He and each valence neutron. The interactions V an 
and V nn are given by the microscopic KKNN potential and the effective Minnesota 
potential, respectively. These potentials well reproduce the scattering data of 4 He-n 
and n-n systems. In this three-body model, there is a small deficiency of the binding 
energy (~ a few hundred keV) of 6 He ground state, which is considered to come from 
the 4 He core polarization effect. 17 ) In order to improve this deficiency of the binding 
energy, we employ the effective ann three-body interaction V^ nn as 

VLn = V,e-^ + ^\ (2-3) 

where V 3 = -1.503 MeV and v = 0.07/6;? far 2 . 

The component of the Pauli forbidden state is excluded in the relative wave 
function x J * \nn) by using the so-called pseudo potential Vp-p. In the case of 6 He 
with the 4 He core, the Pauli forbidden state c/>pf for the valence neutrons is the 
occupied 0s state of 4 He. Then, the pseudo potential Vpp is given as 

2 

F PF = ^A|0 l PF )(4 F |, (2-4) 
i=i 

where A is taken as 10 6 MeV and i is an index for valence neutrons. 

Equation f)2 • 1 j) is solved accurately in a few-body technique. We here employ the 
variational hybrid- TV model, in which the relative wave function of the 4 He+n+n 
system are expanded on the superposed basis states of the cluster orbital shell model 
(COSM; V-basis) and the extended cluster model (ECM; T-basis); 17),22) 

X JV (nn) = X f (£y) + xF (£r), (2-5) 

where Xvt(£vt) expresses the relative wave function, and £ v and £ T are V- and 
T-type coordinate sets, respectively. The radial component of each relative wave 
function is expanded by Gaussian basis functions (Gaussian expansion method, 
GEM). 22 )' 23 ) This model successfully describes the observed properties of 6 He such 
as the two-neutron binding energy (0.975 MeV) and the matter radius (2.46 fm) of 
the + ground state. 

2.2. Complex scaling method 

In CSM, 24 ) -26 ) relative coordinates for a many-body system are transformed as 

U(6)TU- 1 (O) = re i0 , (2-6) 

where U(6) is a complex scaling operator and 9 is a scaling angle given in a real num- 
ber. Applying this transformation to the Hamiltonian H, we obtain the complex- 
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scaled Hamiltonian H e . For H e , the corresponding complex-scaled Schrodinger equa- 
tion is expressed as 

H°X e = E X e , X 9 = e^ e -f X (re* e ), (2-7) 

where x d is a complex-scaled wave function. The factor e^ 3 / 2 ^ 6 'f comes from a 
Jacobian for a volume integral with f degrees of freedom of a system (/ = 2 for a 
three-body system). 

We obtain eigenstates (their biorthogonal states) and energy eigenvalues of the 
complex-scaled Hamiltonian H as {x n } ({XnD an d {E^} with a state index n, re- 
spectively, by solving the eigenvalue problem of Eq. (|2-7|) using a finite number of I? 
basis functions. In CSM, all energy eigenvalues of unbound states are obtained on the 
lower half of a complex energy plane, governed by ABC-theorem, 24 ) and their imag- 
inary parts represent outgoing boundary conditions. In ABC-theorem, it is proved 
that a divergent behavior at an asymptotic region of resonances is transformed to a 
dumping one by CSM. This condition enables us to obtain many-body resonances 
by the same calculational way as the bound state case. The resonances are obtained 
with the complex energy eigenvalues of E d n = E r n — ir n /2, where E r n and F n are 
resonance energies measured from the threshold and decay widths, respectively, and 
these energy eigenvalues are independent of the scaling angle 9. On the contrary, 
energy eigenvalues of continuum states are obtained on the branch cuts of the Rie- 
mann sheet, which are rotated down by 29. This difference of the behaviors between 
resonances and continuum states makes the energy eigenvalues of resonances isolated 
from continuum states as shown in Fig. [TJ Actually, CSM has been widely employed 
as a useful tool to identify the resonance of two-, three- and four-body systems. 3 ^' 27 ^ 

Moreover, CSM is also useful to solve a decay problem of a many-body system. In 
CSM, the obtained energy eigenvalues of continuum states in a three-body system 
are classified into two- and three-body ones by ABC-theorem. These continuum 
states are located on the 20-rotated branch cuts starting from different thresholds 
of two- and three-body decay channels, such as 5 He+n and 4 He+n+n in the case 
of 6 He. (See Fig. [TJ) The classification of continuum states in CSM imposes that 
an outgoing boundary condition for each open channel is taken into account by 
imaginary parts of energy eigenvalues. Using this classification of continuum states, 
we can describe three-body scattering states without any explicit enforcement of 
boundary conditions. 

2.3. Complex-scaled solutions of the Lippmann-Schwinger equation 

We explain a new method of CSLS to describe three-body scattering states, 
which is capable of calculating physical quantities as functions of subsystem energies 
in a three-body system. 

The formal solution of the Lippmann-Schwinger equation can be described as 

tp-W = <p + lim 1 V<P , (2-8) 

^E-H±ie 

where is a solution of an asymptotic Hamiltonian Hq. The total Hamiltonian H 
is given in Eq. (|2-2p for 6 He and the interaction V is given by subtracting Hq from 
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Re(E) [MeV] 

Fig. 1. Energy spectra of 6 He in CSM. The left panel is a schematic picture of spectra of the 
4 He+n+n system, where the indices Ci, C2 and C3 indicate the 4 He+n+n, 5 He(3/2 _ )+7i and 
5 He(l/2~)+n continuum spectra, respectively. The right panel is the 1~ spectra calculated 
within the present 4 He+n+ra model. Three straight lines show the three- and two-body con- 
tinuum states of 4 He+n+n, 5 He(3/2~)+n and 5 He(l/2~)+n corresponding to Ci, C2 and C3, 
respectively. In the 1~ spectra calculated with the scaling angle 8 — 27 degrees, no bound and 
resonant states are obtained. 



H. The boundary condition of the scattering state is represented by ±e. 

When we consider the scattering states of the Borromean system, the asymptotic 
Hamiltonian Hq is equivalent to a free Hamiltonian in a three-body system since all 
the scattering states are described by three-body scattering states and any binary 
subsystem does not have a bound state. Then, we obtain the following relations: 

3 

£ = 5>-2b (2-9) 

i=i 

<?>o = ^3yexp(zk-r)exp(iK-R), (2-10) 

where k, K and r, R are relative momenta and relative coordinates in a three- 
body Jacobi coordinate system, respectively. We denote ^0 by (k,K| and |k, K) 
in the bra- and ket-state representation, respectively, to describe the momenta in 
asymptotic region, k and K, explicitly 

The formal solution in Eq. (|2-8p for an outgoing three-body scattering state can 
be rewritten in the ket-state representation as 

|^ (+) (k,K)) = |k,K) + lim 1 V\k,K), (2-11) 

E - H + ie 

where the interaction V is 

2 

V = H - H = V ^i) + Vnn + VLn + ^PF- (242) 
i=l 
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In the present work, we use the complex-scaled Green's function Q £'). 
The complex-scaled Green's function is related to the non-scaled Green's function 
G(E;t,£') as 

iim - I , = O = u-\e)g 6 (E;t,t>)u(e), (2-13) 

E - H + ie 

where the complex-scaled Green's function G e (E;£,£') is defined as 

1 



E — H e 



j-i \ \ r Xn (flj&(0 



In derivation of the right hand side of Eq. (|2-14|) . we use the extended complete- 
ness relation (ECR), whose detailed explanation is given in Ref. 28) and skipped 
here. Using this Green's function constructed by the complex-scaled wave functions 
of bound, resonant and non-resonant continuum states, we can take into account 
boundary conditions for all open channels of a three-body system in the form of 
complex energy eigenvalues E® n . Then, we can omit the operation of e — > in the 
derivation of the complex-scaled Green's function in Eq. f)2- 14j) . It should be noted 
that the Green's function in Eq. (|2-14p is a continuous function with respect to the 
total energy E while we use discretized energy eigenvalues E e n of the complex-scaled 
Hamiltonian H e . In this calculation, we use 15 Gaussian bases for each relative 
coordinate, the range of which is taken up to about 20 fm, and the scaling angle 9 
is taken as 18 degrees to construct the Green's function. 

Combined with the complex-scaled Green's function in Eq. (|2- 14|) . the formal 
solution in Eq. (|2-11|) is rewritten as 

|^+)(k,K)) = |k,K) + ^C/- 1 (0)|x^) ¥ ^(^l^)l>|k,K}. (2-15) 

It is not necessary to apply the complex scaling to the first term of the solution 
of an asymptotic Hamiltonian. The complex scaling operator U(8) in Eq. (|2-15|) 
is processed in the calculation of the matrix elements and does not operate on the 
complex-scaled eigenstates \n- 

Similarly, let us consider the formal solution for an incoming scattering state. 
To describe the incoming scattering state, we here consider the bra-state of 
with momenta (k, K), which is given as 



(# (_) (k,K)| =(k,K| ( 1 + lim V | 

V e ~*° E-H -ie ) 



= (k,K| + (k,K|T>lim ' 



t 

(2-16) 



In the derivation of the second line, we assume the hermiticity of H and V. The 
Green's function in Eq. (|2- 16j) is equal to that of Eq. (|2-13p . and hence we replce this 
Green's function in to the complex-scaled Green's function. Using Eqs. (|2- 13[) and 
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(I2-I4j) . we obtain the bra-state of the incoming scattering state as 

(ffH(k,K)| = (k,K| + fcK\VU-\0)\£) w ^(&\U(p). (2-17) 

Hereafter, we refer Eqs. (|2-15p . (]2- 17H and their conjugate states as the complex- 
scaled solutions of Lippmann-Schwinger equation (CSLS). 

Additionally, we switch off the pseudo potential Vpp in the calculation of Eqs. f|2- 15[) 
and (|2-17p to avoid an instability of numerical results, which comes from the large 
value of A = 10 6 MeV, while the wave functions x n are solved with the pseudo poten- 
tial. By switching off the pseudo potential, the antisymmetrization in the scattering 
state is approximately ignored, but it is not a serious problem, since the antisym- 
metrization in the intermediate states are considered when we solve the eigenstates 
Xn- In fact, it will be shown in the next section that the obtained result in CSLS 
shows a reasonable agreement with the previous result in Ref. 17) and the calculated 
breakup cross section well reproduces the trend of the experimental data. 

2.4. El transition in complex- scaled solutions of the Lippmann-Schwinger equation 

To calculate the El transition strength and the Coulomb breakup cross sec- 
tion, we start with the two-dimensional momentum distribution of the EX transition 
strength given as 

where <£ g . s . is a ground state wave function and O(EX) is a transition operator with 
a rank A. J g . s . is a total spin of the ground state. 

Using Eq. (|2-18p . we derive an EX transition strength with respect to the total 
energy E of a system as follows; 

dB(E\) _ f [ dMK fB(EX) s ( E _WW\ 



dE JJ dkdK \ 2/2 2M 

where \x and M are reduced masses of subsystems corresponding to the two mo- 
menta k and K, respectively. Similarly, we obtained the two-dimensional energy 
distribution as 

*»™ - it^fmm, u _ **\ , l - ^) , (2.20) 



deide 2 J J dkdK \ 2^ J \ 2M 

where £\ and £2 are subsystem energies in a three-body system. 

In the next section, we shall show the total energy and two-dimensional energy 
distributions of the El transition strength, and discuss the correlations of subsystems 
in the Coulomb breakup of 6 He. 

§3. Energy distributions of El transition strength for 6 He 



We demonstrate that CSLS is useful to investigate internal correlations of sub- 
systems in the final states. 
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Energy [MeV] 

Fig. 3. (Color online) The obtained Coulomb breakup cross section measured from the 4 He+ra+n 
threshold energy. The red and blue lines show the original result in CSLS and the folded one 
with the experimental resolution, respectively. The observed cross section (open square) is taken 
from Ref. 11). 

Before discussing the correlations of subsystems in the final states, we first cal- 
culate the total energy distributions of the El transition strength from the ground 
state of 6 He into 4 He+n+n three-body scattering states to show the importance of 
the final state interaction (FSI). Using Eq. (|2-19p . we obtain the result of the to- 
tal energy distribution measured from the 4 He-|-n-|-n threshold energy as shown in 
Fig. [2j From this original result, it is confirmed that the strength possesses a sharp 
peak at around 1 MeV. In order to recognize FSI, we also calculate the strength 
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without FSI. When we switch off FSI, scattering state of He can be described only 
by first term of the right hand side of Eq. (|2-15p since the interaction V is zero. Then, 
the strength without FSI, which is equivalent to the one of the transition from the 
ground state into non-interacting three-body continuum states, is calculated by tak- 
ing the first term of Eq. ()2- 15|) . This transition strength of the direct breakup is also 
shown in Fig. [2l It is found that the direct breakup one has a small strength with 
a broad bump at around 3 MeV. This large difference between the original and the 
direct breakup strengths indicates the importance of FSI to explain the low energy 
enhancement in the El transition strength of 6 He. The result of Fig. [2] is consistent 
with the previous result in Ref. 17). 

Here, we also show the reliability of our calculation by comparing the obtained 
result using CSLS and the experimental data. We drive the breakup cross section of 
6 He using the obtained El transition strength and the equivalent photon method. 
The experimental data are taken from Ref. 11). In Fig. [3l two strengths are shown; 
One is the original result in CSLS and another is the folded one by the experimental 
resolution. 11 ^ The obtained cross section in CSLS has a peak at around 1 MeV and 
well reproduces the observed trend. This good agreement of the cross section implies 
the reliability of CSLS to describe the three-body scattering states since the validity 
of the three-body model has been already shown in Ref. 17). 

Next, we investigate correlations of subsystems in the final states. Using Eq. (|2-20p . 
we evaluate the two-dimensional energy distribution of the El transition strength of 
6 He, associated with the 5 He subsystem. The result is shown in Fig. [H where E\ and 
£2 in Eq. (|2-20p are the relative energies of 4 He-n (E a . n ) and 5 He-n (Es^ e _ n ) systems, 
respectively. It is clearly seen that the strength is concentrated on around E a _ n ~ 0.7 
MeV, which agrees with the 5 He(3/2~) resonance energy. Hence, the importance of 
the 5 He(3/2 _ ) resonance in the Coulomb breakup reaction of 6 He is directly shown 
in the physical observables using CSLS, which is consistent with the previous anal- 
ysis using RFM. 17 ) We also show the two-dimensional energy distribution of the 
direct breakup strength in Fig. [5] to clarify the FSI in the two-dimensional energy 
distributions. From Fig.[5l we find that no clear peak structure appears without FSI. 
This result indicates that the correlations in the two-dimensional energy distribution 
mainly comes from the FSI and the sequential decay via the 5 He(3/2 _ )+n channel 
is important in the Coulomb breakup reaction of 6 He. 

§4. Summary 

In the present study, we developed a new approach called the complex-scaled 
solutions of Lippmann-Schwinger equation (CSLS), which enables us to describe a 
scattering state for three-body breakups of the Borromean system and to calculate 
the observables with respect to not only the total energy but also the subsystem 
energies. Using CSLS, we reproduced the observed Coulomb breakup cross section 
nicely, and confirmed that the 5 He(3/2~) resonance plays an important role in the 
Coulomb breakup reaction. This CSLS approach makes us to extract the correlations 
of subsystems from the observables and is useful to study the properties of the 
weakly-bound nuclei not only for the ground state, but also for the scattering states. 
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0.6 1.2 1.8 2.4 3 
E a -n [MeV] 

Fig. 4. (Color online) Twi-dimensional energy distribution of the El transition strength corre- 
sponding to the 4 He-n and 5 He-n subsystems. 




0.6 1.2 1.8 2.4 3 
E a -n [MeV] 

Fig. 5. (color online) Two-dimensional energy distribution of the direct breakup as same as Fig. [4] 



The detailed analysis on the subsystems correlations such as He-n and n-n systems 
in the Coulomb breakup reaction of 6 He is forthcoming. It is also interesting to 
perform the analysis of the two- neutron halo nuclei 11 Li using this method. 
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